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Abstract 

We discuss distortion in microlensing-induced light curves which are considered to 
be curves due to single-point-mass lenses at a first glance. As factors of the distortion, 
we consider close binary and planetary systems, which are the limiting cases of two- 
point-mass lenses, and the gravitational potential is regarded as the sum of the single 
point mass and the corrections. In order to quantitatively estimate the asymmetric 
features of such distorted light curves, we propose a cutoff dependent skewness, and 
show how to discriminate the similar light curves with it. We also examine as the 
distortion the general relativistic effect of frame dragging, but the effect coincides with 
the close binary case in the light curves. 
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§1. Introduction 

Gravitational microlensingi^ is one important probe for studying the nature and distri- 
bution of mass in the galaxy. Chang and RefsdalB and Gottfl suggested that even though 
multiple images by lensing are unresolved, the time variation of the magnitude of the source 
can still be detected if the lens moves relative to the source. Such light curves caused by 
the microlensing can be distinguished from curves of intrinsically variable sources, because 
the change of magnitude by lensing is achromatic, whereas the colors of intrinsically variable 
stars change in general. Furthermore, microlensing-induced light curves can be distinguished 
from magnification by bursts, which are likely to appear in sheer shapes. As is well known, 
microlensing by a point mass has the time-symmetric light curves, provided that the lens 
and the source have constant relative transverse velocity. Hence, most events are expected 
to have almost time-symmetric light curves. If the relative velocity is not constant, then of 
course time-asymmetric light curves will be detected. Gould ^ predicted a parallax effect 
due to the orbital motion of the Earth, and its effect was indeed detected by Alcock el al.E* 
When the time scale of a microlensing event is larger than ~ 100 days, this effect is impor- 
tant. However, if the time scale is of hours to weeks, the parallax effect can be neglected. 
Then, other factors, such as the non-spherical gravitational potential of the lens, could also 
produce distortion from the time-symmetric light curves. The effect may be regarded as a 
higher-order correction to the point-mass lens. The aim of this paper is to evaluate the effect 
of the intrinsic nature of the lenses, which may slightly distort the light curves. We exclude 
significantly peculiar light curves, such as double peaks, from which direct information is 
available without any detailed analyses. We rather restrict our consideration to light curves 
which are regarded as curves of single-point-mass lenses at a first glance, that is, almost time- 
symmetric light curves. Since the asymmetric part contains additional information about 
massive astrophysical compact halo objects (MACHOs), this subject is very important to 
understanding the nature of MACHOs. 

One of the important factors to induce time-asymmetric forms is the binary system of the 
lensing objects, in which the contribution from both objects to the Newtonian gravitational 
potential is no longer spherical. The discussion is simplified by considering two-point-mass 
lenses. In particular, since we consider almost time- symmetric light curves, we hit on two 
situations: close binary and planetary systems. In a close binary system, the separation 
distance between the two point masses is much less than the Einstein radius of the total 
mass, so that the approximated light curve can be described by the total mass and some 
corrections to it. The detectability of a close binary system has been discussed in detail by 
Gaudi and Gould0> considering the excess magnification threshold. According to them, the 
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detect ability of close binaries with separation less than ~ 0.2 of the Einstein radius is ~ 10%. 
Therefore, most close binary lenses with such small separation are missed. We consider the 
possibility of picking up the discarded asymmetry, by which even below the separation of 
~ 0.2 of the Einstein radius the binary nature of the lenses may be detected. In a planetary 
system, the light curve is expected to be described by the contribution from the larger mass 
and some corrections from the smaller one. Several authors0>J)>i)'0><0> have discussed plan- 
etary systems with remarkable deviations from a single lens. However, we are now interested 
in planetary events which would be missed due to the special configurations of the lensing 
geometry. Our new proposal may make the detection of such events possible, in addition to 
the close binary case. Another factor for the time-asymmetric forms is the asymmetry of the 
lens object itself. The multipole moment of the lens deviates from gravitational potential 
of the point mass. The gravitational potential due to the quadrupole never vanishes for the 
two-point-mass case. Therefore, such a correction term in the gravitational potential can 
be regarded as the limiting case of the binary. Instead, we shall consider general relativistic 
effects of dragging of inertial frames due to a rotating object as another factor, which cannot 
be expressed as corrections to the Newton-like potential. While the Newton-like potential 
corresponds to the gravitoelectric field, this effect results from the gravitomagnetic field. 
In this case, corrections up to post-Newtonian order are sufficient. Other post-Newtonian 
potentials also affect the light curves, but they are neglected here. 

This paper is organized as follows. In §Q we discuss the two situations involving the two- 
point-mass lenses and obtain microlensing- induced light curves slightly deviating from the 
time-symmetric curves. In §|3] we discuss the post-Newtonian corrections due to the rotation 
and obtain time-asymmetric light curves as well. We propose a certain tool to estimate the 
time-asymmetric features quantitatively in §[|. Finally, we give summary and discussion in 
§1- 



§2. Two-point-mass lenses 

We investigate two extreme situations involving two-point-mass lenses in order. (The 
general computational details including critical lines and caustics can be found in Ref. |12|).) 
First, we consider the situation in which the distance I between the two point masses is 
much less than the Einstein radius te corresponding to the total mass. Next, we consider 
the situation in which one mass, Mi, is much smaller than the other mass, Mi. 
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2.1. Close binary system 

We consider the lens plane on to which the positions of the two point masses Mi and 
M 2 are projected (see Fig. |l|). We introduce an angular coordinate system (8 X , 9 y ), in which 
the two point masses lie on the # x -axis and the origin is chosen as their geometrical center. 
We also define an angular coordinate system {(3 x ,(3 y ) in the source plane, corresponding to 
the lens plane. We express the distance between the observer and the lens plane, the lens 
plane and the source plane, and the observer and the source plane by Dl, Dls and D$, 
respectively. 

If we denote the angular separation between the source and the image by a, we have 
AGMi D LS e x - V , AGM 2 D LS 6 X + V 
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where t] is the angular separation of the mass Mi (or M 2 ) from the optical axis. Therefore, 
the lens equation can be written in the form 



13 = 6- ex. 

We normalize this equation by the Einstein radius for the total mass, 

(AG {Mi + M 2 ) D LS \ * 



9, 



D L D S , 



(2-2) 



(2-3) 




lens plane 



source plane 



Fig. 1. Geometry of the gravitational lensing considered in §2.1. 
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and introduce normalized quantities 






The lens equation is then given by 



~ @x — V d x + fj 

Px = x -fl 1 — -2 --M27Z T2 -> ( 2 ' 5a ) 

(o,-v) +ei (e x + fj) +ei 

9 6 

Pv = dy-m— v - 2 — --^2— y - 2 — -» ( 2 - 5b ) 

where /xi and /x 2 are defined as 

Mi M 2 
/il " M 1 + M 2 ' ^ " M 1 + M 2 - (2 ' 6) 

The separation distance in the projected plane is I = ■ 2rj, and the Einstein radius for 
this system te = Dl ■ &e- The term 'close binary' in the gravitational lens means I <C te- 
This condition in the astronomical situation is expressed as I <C 1O 14 (M/M ) 1 / 2 (D/(1Q 
kpc)) 1 / 2 cm, where we have chosen typical astronomical distances as the scale of the galactic 
halo, Dls ~ Dl ~ Dg ~ D. Therefore, the range of applicability is not so severely limited. 
The condition of the close binary, I <C Te, is mathematically expressed as 

fj < 1. (2-7) 

Under this condition, we expand the right-hand side of the lens equation with respect to 
fj. Up to first order in fj, we have 

n Q~2 _ q2 

/3 x = 0x- ~ 9 - ' -„ ~ V (pi ~ M2) —r — z\z, (2-8a) 



e l + e l (ei + H 

h = h - Ma) ,- 29x9 H,2 , (2-»b) 



where the last term on each right-hand side represents the deviation from a single-point-mass 
lens. The first-order corrections vanish for the equal mass case, since /ii = /i 2 . It is therefore 
convenient to express the combination fj (/ii — /i 2 ) as one small parameter e = fj (/ix — A^)- 

The inversion of Eq. (|2-8|), i.e., solving 6 by /3, is possible, but the general form is quite 
messy. However, under the condition e C 1, the approximate solution, i.e., the first-order 
solution in e, is given by the form 

= O + (2-9) 
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where 6q and 0± denote the zeroth-order and the first-order solutions, respectively. Substi- 
tuting Eq. ( |2-9| ) into Eqs. ( |2-8| ), we can find such solutions. The zeroth-order solution is 
given by 



®0x — - ( Px =•= Px 



1 + - 



\ Pl + P, 



\ P 2 X + P 



^On 17 - — I Py i Py A 

Using this zeroth-order solution, the first-order solution is written in the form 



AqL ~ °0y ~ 1 



e 



- 1 



(2-10a) 
(2-10b) 

(2-lla) 
(2-llb) 



Next, we turn our attention to the derivation of magnification. The magnification M is 
given by 



M M. ■ M 

'dPi 



det 



09 



3/ + 



+ 



det 



(2-12) 



where the subscript (+) and (— ) correspond to solutions with a plus sign and with a minus 



first order in e, in the following way: 
i -l 



sign, respectively, in Eqs. ( pT0|) . The inverse of the Jacobian det ( if 1 ) is calculated, up to 
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(2-13) 



Therefore, using Eqs. ( |2-1U| ), (2-12) and (2-13), we can derive the magnification as a function 
of the source position: 

M = M(p x ,P y ;e). (2-14) 

The time variation is produced by the relative change of the positions. The time variation 
of the magnitude Am is then given by 

Am = Am{t) = 2.5 log 10 M e) , (2-15) 
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Fig. 2. Microlensing- induced light curves with (solid) and without (dashed) correction for the 
deviation from a single-point-mass lens. The light curve with the correction is plotted for the 
small parameter value e = 0.1. The relative motion of the source to the lens is assumed to be 
described by (p = and p = 0.3. 



where the relative source trajectory is in general described by 

(3 x (t) = — cos (p + p simp, (2- 16a) 

to 

P y (t) = — simp — pcosip, (2-16b) 

where i is the time taken to cross the Einstein radius, tp is the angle of the trajectory from 
the /Jj-axis, and p is the impact parameter normalized by the Einstein radius. From this, we 
can derive the microlensing- induced light curves, slightly deviating from the time-symmetric 
curves. An example of the light curves is shown in Fig. 2 for (p = and p = 0.3. The solid 
line denotes the light curve with the correction for the deviation from a single-point-mass 
lens, while the dashed line corresponds to the single-point-mass lens. The more general 
dependence on the angle <p and the impact parameter p is discussed in detail in §f|. 

2.2. Planetary system 

In this subsection, we discuss the case that the mass Mi is much smaller than M 2 ; that 
is, the object with smaller mass is regarded as the planet. In this situation, it is convenient 
to chose the origin of the (9 X , 8 y ) system to be at the position of mass M 2 . Furthermore, let 
the angular separation between the mass Mi and M 2 be r\ (see Fig. 3). Then the normalized 
lens equation is given by 

k = h - Jj-~- - * _ dx ~J _ , (2.17a) 
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Fig. 3. The lens plane considered in §2.2. The angular separation of the two point masses is 
denoted by r\. 



Py = y 



0>, 



0,, 



>17b) 



where fi = M1/M2 (<C 1), and we have used for the normalization the Einstein radius of the 
mass M2, 

_ (AGM 2 D LS y 

9e -{—d^dS) ■ (248) 

As in the previous subsection, in the situation jii < 1, we approximate the solutions of the 
lens equation up to first order in \i. Using the zeroth-order solution ( |2T0| ), we obtain the 
first-order solution 
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Furthermore, we can calculate the magnification by using Eq. (2-12) and the inverse of 



the Jacobian det ( If 1 ) , which is up to first order in /i, given by 
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(2-20) 



Therefore, we can also obtain the microlensing-induced light curves 

Am = Am(t) =2.5 \og l0 M ((3 x (t) , /3 y (t); fi) , (2-21) 

where the relative source trajectory is described by Eq. ( |2-16| ) using the angle ip and the 
impact parameter p. Some examples of the light curves with different values of fj are shown 
in Fig. 4 for p, = 0.05, ip = and p = 0.3. As seen in these figures, some light curves tend 
to have double peaks with certain geometrical configurations. The second peak corresponds 
to the effect of a planet of mass 0.05M2. The ratio of the height at the peaks is almost 
determined by the mass ratio fi. Although we have restricted ourselves to the case /i < 1 
in order to exclude peculiar light curves, light curves with double peaks are still obtained 



'0x 



v) + Oly 



becomes 



with certain configurations. This is because the correction term 
effective when 9 Qx is equal to fj. (More detailed discussions of planetary-binary lensing with 
dramatic features are given in Refs. f7|)-|ll]).) 

Furthermore, it is interesting to investigate the configurations in which the Lorentzian 
curves, i.e., almost time-symmetric light curves arise. For this purpose, we consider the 
integral of the square of the magnitude difference from the corresponding single-point-mass 
lens: S = f_™(Am — Amo) 2 dt. The quantity S indicates the criterion of the deviation from 
the light curve due to the single-point- mass lens. Figure 5 displays the contours of 5 in the 
fj-tp space for different impact parameters. Of course, as the quantity S becomes smaller, 
the light curve moves closer to that due to the single-pint-mass lens. In fact, we can find 
almost time-symmetric light curves in the domains where 5 is smaller than ~ 0.005, as seen 
in Fig. 4(b). Such domains tend to become larger as the impact parameter p increases. The 
same tendency can also be derived by making the mass ratio \x small. 



§3. Rotating objects 

Several relativistic effects also causes corrections to the point-mass lens. We only consider 
the dragging effect of inertial frames arising from a rotating object, since other spherical post- 
Newtonian terms never induce asymmetry in light curves. This additional effect is described 
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Fig. 4. Microlensing-induced light curves with (solid) and without (dashed) correction for the 
deviation from a single-point-mass lens. The light curves with the correction are plotted for 
the case of the small parameter value = 0.05 and for angular separations of (a) fj = 0.3, 
(b) fj = 1.0, and (c) i) = 1.7. The relative motion of the source to the lens is assumed to be 
described by ip = and p = 0.3. 
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. 5. Contours of 5 in 77-^ space for the impact parameters p = 0.2, 0.3 and 0.4. The small 
parameter ji is assumed to be 0.05. The attached labels 'a'-'r' indicate the parameters used in 
Fig. 10 (see text). 
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Fig. 6. Geometry of lensing by the rotating object considered in §3. 



with the spin angular momentum J. We consider the projection of the angular momentum 
of the rotating object on to the lens plane and define the (6 X , 9 y ) coordinate system so that 
the ^-axis is oriented parallel to the projected angular momentum J± (see Fig. |]). The 
deflection angles (a x ,a y ) are written as© 

agm e x agj ± e 2 x -e 2 y 

ax ~ Jn Q2 , Q2 + Jn2 ~T~ ~V2' (3-la) 



c 2 D L e 2 x + e 2 y c 3 d 2 l 



Hence, we have 



AGM 9 y AGJ ± 26Jy 

ay -cW- L WTo 2 + c*d 2 ^ x + el y (34b) 
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Therefore, using the quantities normalized by the Einstein radius, the lens equation becomes 
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4 = ^-^-77^^, (3-3b) 
°* + °l {01 + ft, 2 ) 

where 7 is given by 

It is interesting that this lens equation is the same as that of the two-point-mass lenses for the 
close binary system, I <C te, with the correspondence 7 <-> s. Hence, the same asymmetric 
light curves are obtained. Furthermore, it is impossible to distinguish two such corrections. 
However, the parameter 7 is quite small. For example, we consider a Kerr black hole. The 
angular momentum is J ~ GM 2 /c, so that we have 7 ~ (GM D s j \c 2 D L D ls)) 1 ^ 2 "C 1. 



§4. Quantitative estimate 

In order to estimate the asymmetry in the light curves quantitatively, we now introduce 
the notion of 'skewness' from statistics. 0* In statistics, the skewness for any distribution 
function f(t) is defined as 



skewness 



1 

N 



00 ft - yU 



a 



f(t)dt, 



(4-1) 



where N is the normalization factor, \x the mean, and a the standard deviation. However, 
there is a problem in utilizing the skewness exactly in this form. To see this, we consider a 
single-point-mass lens. The light curve is then given by 



Am = 2.5 log 
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Fig. 7. Illustration of the discussion given in §4. A schematic light curve with noise. 
From this, we find that the integral 



observed light curves necessarily include noise, it seems that the integral to infinity is mean- 
ingless. With this consideration, we define the skewness for the restricted region of a light 
curve Am > AZ\mmax , where Amm&x is the maximum value. Here we have introduced the 
cutoff A (0 < A < 1), and the usual skewness corresponds to A = 0. The cutoff will naturally 
appear in the actual data, e.g., the region Am < AZ\m max may be meaningless due to the 
noise level (see Fig. 0). In microlensing-induced light curves, the bottom level is constant, 
but the maximum value may have some uncertainty as (1 ± A)Z\mmax- 

We apply this tool to the case of almost identical light curves, which are caused from quite 
different physical situations, that is, a close binary and a planetary system. The skewness 
of these light curves are respectively displayed as functions of 1/A in Figs. 8-11. Figures 8 
and 9 display the dependence of the skewness of the close binary light curves on the angle 
ip and the impact parameter p, respectively. The light curve given in Fig. 2 corresponds to 
the curves in Figs. 8 and 9. On the other hand, Figs. 10 and 11 display the skewness of the 
light curves in the planetary system. In the case of the planetary system, the almost time- 
symmetric light curves are derived only under certain configurations, as indicated in §2.2. 
Figure 10 displays the skewness of the light curves under such configurations. The curves 
labeled by 'a'-'r' correspond to the points labeled by 'a'-'r', respectively, in Fig. 5. Figure 
11 displays the impact parameter dependence of the curves corresponding to the point 'a'. 
The light curve given in Fig. 4(b) corresponds to the curve 'a' in Fig. 10 and a curve in 
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diverges if n > 3. Therefore, we cannot 
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Fig. 11. 
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Fig. 8. Dependence of the skewness on the angle <p in the close binary case. The trajectories have 
the same impact parameter, p = 0.3, and the light curves are calculated for the small parameter 
value e = 0.1. 



As is expected, we have zero skewness when (p = ±| in Fig. 8. However, when tp is 
different from ±|, we have comparable values of the skewness. This demonstrates the useful 
aspect of the method using the skewness. Nevertheless, the skewness becomes smaller as 
the impact parameter increases, as seen in Figs. 9 and 11. Thus, the usefulness of the 
method using the skewness depends mainly on the impact parameter. The skewness with 
respect to different values of A fully shows the asymmetric features of the light curves. The 
absolute values of the skewness depend on the the small parameters e and /x (i.e., the binary 
separation and the mass ratio). Furthermore, it should be noted that the absolute values 
of the skewness have a maximum at 1/A ~ 2-4 and decrease monotonously in the close 
binary case, while in the planetary case the skewness indicates different behavior, peaks at 
larger values of 1/A, and so on. Therefore, we may discriminate the underlying cases for the 
asymmetry by the A-dependent skewness for a good signal-to-noise ratio beyond ~ 10. 



§5. Summary and discussion 

We have studied distortion in microlensing-induced light curves which are considered to 
be curves of single-point-mass lenses at a first glance. In particular, we have attributed 
factors inducing the distortion to lenses themselves and considered two sorts of corrections: 
corrections due to deviation from the Newtonian gravitational potential = —GM/r and 
corrections due to general relativistic effects of dragging of inertial frames arising from a 
rotating object. For simplicity, we have discussed two extreme cases in two-point-mass 
lenses for the corrections of the potential; one is the close binary case in which I «r B , and 
the other is the planetary system case in which M 1 <C M 2 . Moreover, we have considered 
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Fig. 9. Dependence of the skewness on the impact parameter p in the close binary case. The 
trajectories have the same angle, ip = 0, and the light curves are calculated for the small 
parameter value e = 0.1. 

corrections up to the post-Newtonian order for the effect of dragging of inertial frames. From 
this, we found the same time-asymmetric light curves as in the two-point-mass lenses where 
I <C r^. Furthermore, we have introduced the cutoff dependent skewness and estimated the 
asymmetry in the light curves quantitatively. In particular, we showed the clear difference 
in the skewness for almost similar light curves. 

Here we make a comment on the additional factor for the asymmetry in the binary 
system. In this paper, we have assumed that the lens systems of the two point masses are 
fixed, but each star constituting the binary revolves around the center of mass. Therefore, the 
rotational effect may also appear. However, if the rotation period T ~ (l 3 /(G(M 1 + M 2 )) 1 ^ 2 
of the binary is much larger than the typical time scale to of a microlensing event, our 
consideration of the fixed lens systems would be appropriate. For the extreme case T <^t , 
the time-averaged gravitational potential which is projected on to the lens plane may be 
regarded approximately as that of a single-point-mass lens or a fixed two-point-mass lens 
with I <C r E if the lens is compact. Therefore, more complicated variations of light curves, 
corresponding to the phase, are expected only if T ~ t . 

It is very interesting whether the time-asymmetric features of light curves discussed in 
this paper will actually be detected by the projects (MOA, etc.) that are now in progress. 
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Fig. 10. The skewness of light curves in the planetary system. The curves labeled 'a'-'r' correspond 
to the points labeled 'a'-'r', respectively, in Fig. 5. The light curves are calculated for the small 
parameter value fi = 0.05. 



References 

[1] See, e.g., B. Paczyrisky, Ann. Rev. Astr. Ap. 34 (1996), 419; R. Narayan and 

M. Bartelmann, Lectures on Gravitational Lensing, preprint: |astro-ph / 960600 1| . 
[2] K. Chang and S. Refsdal, Nature 282 (1979), 561. 
[3] J. R. Gott, Astrophys. J. 243 (1981), 140. 



17 



-0.02 1 ■ ■ — ■ ■ — ■ J 

1 10 100 

1 IX 

Fig. 11. Dependence of the skewness on the impact parameter p in the planetary case. The light 
curves correspond to the point 'a' in Fig. 5, and are calculated for the small parameter value 
H = 0.05. 

[4] A. Gould, Astrophys. J. 392 (1992), 442. 

[5] C. Alcock, R. A. Allsman, D. Alves, et al, Astrophys. J. 454 (1995), L125. 

[6] B. S. Gaudi and A. Gould, Astrophys. J. 482 (1997), 83. 

[7] S. Mao and B. Paczyhski, Astrophys. J. 374 (1991), L37. 

[8] A. Gould and A. Loeb, Astrophys. J. 396 (1992), 104. 

[9] A. D. Bolatto and E. E. Falco, Astrophys. J. 436 (1994), 112. 
[10] D. P. Bennett and S. H. Rhie, Astrophys. J. 472 (1996), 660. 
[11] B. S. Gaudi and A. Gould, Astrophys. J. 486 (1997), 85. 
[12] P. Schneider and A. Weiss, Astron. Astrophys. 164 (1986), 237. 
[13] S. O. Sari, Astrophys. J. 462 (1996), 110. 

[14] See, e.g., W. H. Press et al., Numerical Recipes in Fortran (Cambridge University 
Press, 1992). 



18 



